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Abstract 

We assume that the level spectra of quantum systems in the initial phase of transition from 
integrability to chaos are approximated by superpositions of independent sequences. Each indi- 
vidual sequence is modeled by a random matrix ensemble. We obtain analytical expressions for 
the level spacing distribution and level number variance for such a system. These expressions are 
successfully applied to the analysis of the resonance spectrum in a nearly integrable microwave 
billiard. 
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Random matrix theory provides a framework for describing the statistical properties of spec- 
tra for quantum systems whose classical counterpart is chaotic . It models the Hamil- 
tonian of the system by an ensemble of random matrices, subject to some general symmetry 
constraints. Time-reversal-invariant quantum system are represented by a Gaussian orthog- 
onal ensemble (GOE) of random matrices when the system has rotational symmetry and 
by a Gaussian symplectic ensemble otherwise. Chaotic systems without time reversal sym- 
metry are represented by the Gaussian unitary ensemble (GUE). A complete discussion of 
the level correlations even for these three canonical ensembles is a difficult task. Most of 
the interesting results are obtained for the limit of large matrices. Analytical results have 
long ago been obtained for the case of two-dimensional matrices J3]. It yields simple analyti- 
cal expressions for the nearest-neighbor-spacing distribution (NNSD), renormalized to make 
the mean spacing equal one. The spacing distribution for the GOE, p(s) = |sexp (— fs 2 ), 
where s is the spacing between adjacent energy levels rescaled to unit mean spacing D, is 
known as Wigner's surmise. Analogous expression p(s) = ^§s 2 exp (— ^s 2 ) , is obtained for 
the GUE 

There are elaborate theoretical arguments by Berry and Tabor 5] that classically inte- 
grable systems should have Poissonian statistics. The Poisson distribution of the regular 
spectra has been proved in some cases (see, results by Sinai [s| and Marklof [zj], for in- 
stance). Still its mechanism is not completely understood. It has also been confirmed by 
many numerical studies, although the deviations of the calculated P{s) from exp(— s) are 
often statistically significant (see 8|] and references therein). The appearance of the Pois- 
son distribution is now admitted as a universal phenomenon in generic integrable quantum 
systems. 

A typical Hamiltonian system shows a phase space in which regions of regular motion 
and chaotic dynamics coexist. These systems are known as mixed systems. Their dynamical 
behavior is by no means universal, as is the case for fully regular and fully chaotic systems. 
If we perturb an integrable system, most of the periodic orbits on tori with rational fre- 
quencies disappear. However, some of these orbits persist. Elliptic periodic orbits appear 
surrounded by islands. They correspond to librational motions around these periodic or- 
bits and reflect their stability. The Kolmogorov-Arnold-Moser (KAM) theorem states that 
invariant tori with a sufficiently incommensurate frequency vector are stable with respect 
to small perturbations. Numerical simulations show that when the perturbation increases 
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more and more tori are destroyed. For large enough perturbations, there are locally no tori 
in the considered region of phase space. The break-up of invariant tori leads to a loss of 
stability of the system, that is, to chaos. There are three main scenaria of transition to 
global chaos in finite-dimensional (nonextended) dynamical systems, one via a cascade of 
period-doubling bifurcations, a Lorenz system-like transition via Hopf and Shil'nikov bifur- 
cations, and the transition to chaos via intermittence , 111] . It is natural to expect that 
there are other (presumably many more) such scenaria in extended (infinite-dimensional) 
dynamical systems. 

The nature of the stochastic transition is more obscure in quantum than in classical me- 
chanics. So far in the literature, there is no rigorous statistical description for the transition 
from integrability to chaos. The assumptions that lead to the RMT description do not apply 
to mixed systems. While some elements of the Hamiltonian of a typical mixed system could 
be described as randomly distributed, the others would be non-random. Moreover, the ma- 
trix elements need not all have the same distributions and may or may not be correlated. 
Thus, the RMT approach is a difficult route to follow. Comprehensive semiclassical compu- 
tations have been carried out for Hamiltonian quantum systems, which on the classical level 
have a mixed phase space dynamics (see, e.g. [l2|] and references therein). Berry and Rob- 
nik elaborated a NNSD for mixed systems based on the assumption that semiclassically the 
eigenfunctions and associated Wigner distributions are localized either in classically regular 
or chaotic regions in phase space [13j. Accordingly, the sequences of eigenvalues connected 
with these regions are assumed to be statistically independent, and their mean spacing is 
determined by the invariant measure of the corresponding regions in phase space. There 
have been several proposals for phenomenological random matrix theories that interpolate 
between the Wigner-Dyson RMT and banded random matrices with an almost Poissonian 
spectral statistics [l4]. Unfortunately, these works do not lead to valid analytical results, 
which makes them difficult to use in the analysis of experimental data. There are other 
phenomenological approaches (see, e.g. [ijj] and references therein), which use nonexten- 



sive statistical mechanics, based on maximizing Tsallis or Kaniadakis entropies [16l Il7j , as 
well as the recently proposed concept of superstatistics. These approaches have the advan- 
tage of conserving base invariance of the Hamiltonian matrix. They provide a satisfactory 
description near the end of transition from integrability to chaos. 

This paper considers another phenomenological approach, which has the spirit of the 
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KAM theorem, to the stochastic transition in quantum systems. The phase space of the 
integrable system consists of infinitely many tori corresponding to the conserved symmetries 
of the system. In the semiclassical limit, energy eigenstates are expected to be localized on 
individual tori. Tori destruction corresponds to the mixing of the corresponding quantum 
eigenstates. Symmetry breaking breaks some of the invariant tori but only deforms others 
according to the KAM theorem. Quantum symmetry breaking strongly mixes a limited 
number of eigenstates, but has a less influence on the other ones. Thus, the spectrum is di- 
vided into independent sequences of eigenvalues. States belonging to the same sequence are 
strongly mixed. The sequence may be modelled by a GOE if time reversal invariance is pre- 
served. The interaction between states belonging to different sequences grows as symmetry 
breaking increases. This amounts to amalgamating the initial sequences into a fewer number 
of independent sequences with no more regular character. Consequently, as the number of 
the no-symmetry sequences decreases, their fractional density increases accordingly. As the 
state of chaos is reached, the whole spectrum consists of a single (GOE) sequence. 

Abul-Magd and Simbel 18] consider another class of mixed systems, in which the degrees 
of freedom are divided into two noninteracting groups, one having chaotic dynamics and one 
regular. The Hamiltonian of such a system is given as a sum of two terms, so that each of 
the eigenvalues of the total Hamiltonian is expressed as a superposition of two eigenvalues 
corresponding to the two Hamiltonian terms. The spectrum is then given by a superposition 
of independent chaotic subspectra. Each subspectum corresponds to one (or one set) of the 
quantum numbers of the regular component of the Hamiltonian. This model is used in 
to describe level statistics of vibrational nuclei. An elaborated version of this model 2(| has 
been applied to study NNSD of a wide range of nuclei [211 ]. 

We shall now consider the energy spectra of nearly integrable systems that may be rep- 
resented as a superposition of independent sequences Si each having fractional level density 

m 

fi , with % = 1, ...,m, and with ^ fi = 1. In this case, NNSD of the composite spectrum 

i=i 

can be exactly expressed in terms NNSD's of the constituting sequences (see, e.g., Appendix 
A. 2 of Mehta's book [l|). The gap probability function 

/oo /*oo 
ds' J ds n p{s n ) (1) 

that gives the probability of finding no eigenvalues in segment of length s of the total 
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spectrum, is expressed as a product of the gap functions of the individual sequences 

E(m,s)=l[E i (f i s). (2) 
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We assume that all of Si obey the GOE statistics. Then NNSD of each of the individual 
sequences distribution is given by the Wigner surmise, then for all i 
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where Erfc(x) is the complimentary error function. The NNSD of the full spectrum can the 
be obtained by twice differentiating the resulting gap function. 

We shall restrict our consideration to the case when all sequences have the same fractional 
level density, so that 

fi = f = 1/m. (4) 



The gap function of the composite spectrum is then given by 

E(m, s) 
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Differentiating this function twice with respect to s, we obtain 
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It is easy to see that, for a single sequence p(l,s) = |s exp (— ^s 2 ) , so that the Wigner 
surmise is recovered. On the other hand, lim m ^oop(m, s) = e~ s as required.- 

A weak point of the distribution in Eq. [6] is that it differs from zero at s — 0, because the 
symmetry-breaking interaction lifts the degeneracies. The model thus fails in the domain of 
small spacings as far as the NNSD's are concerned. The magnitude of this domain depends 
on the ratio of the strength of the symmetry-breaking interaction to the mean level spacing. 
Therefore, it is expected to work well for nearly integrable system. 

In the case when the individual sequences are described by a GUE, the gap function for 
each individual sequence is given by 
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In this case, the NNSD of the composite spectrum is given by 

m-2 
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The situation with the level number variance (LNV) E 2 of composite spectra is not as 



clear as in the case of NNSD. Seligman and Verbaarschot [22] argued that S 2 is a variance 
and can therefore be expressed for a composite spectrum as a sum of the corresponding 
quantities for its subspectra, 

m 

S 2 (m,L) = ^S 2 (/ l L), (9) 

where S 2 (x) is the LNV of the ith sequence. There, the LNV of the composite spectrum 
composed of m independent sequences described by RMT is given by 



S 2 (m, L) = m E 2 lMT (L/m) 



(10) 



where S 2 lMT (L) is the LNV calculate by RMT. Explicit expressions for E 2 lMT (L) in the cases 
of GOE and GUE are given in Mehta's book. 

We shall compare our predictions for the NNSD and the LNV with the energy spectra 
of a Limagon billiard. This is a closed billiard whose boundary is defined by the quadratic 
conformal map of the unit circle 2 to w, 



w 



z + Xz , \z\ 



1. 



(11) 



The shape of the billiard is controlled by a single parameter A. For < A < 1/4, the 
Limagon billiard has a continuous and convex boundary with a strictly positive curvature 
and a collection of caustics near the boundary. At A = 1/4, the boundary has zero curvature 
at its point of intersection with the negative real axis, which turns into a discontinuity for 
A > 1/4. The classical dynamics of this system aad the corresponding quantum billiard 
have been extensively investigated by Robnik and collaborators [23]. They concluded that 
the dynamics in the Limagon billiard undergoes a smooth transition from integrable motion 
at A = via a soft chaos KAM regime for < A < 1/4 to a strongly chaotic dynamics for 
A = 1/2. We assume that the quantum dynamics of the Limagon billiard can approximately 



be described by the model present here. The spherical billiard for which A = has two good 
quantum numbers, namely the energy and angular momentum. As A increases, the spheri- 
cal symmetry is gradually destroyed. States corresponding to different angular-momentum 
quantum numbers mix to different degrees depending on the magnitude of the wavefunctions 
at large z. 

The resonance spectra in microwave cavities with the shape of billiards from the family 
of Limagon billiards have been constructed for the values A = 0.125^0.150,0.300 and the 



first 1163, 1173 and 942 eigenvalues were measured, respectively 24|, |25J. The billiard with 



A = 0.300 has a chaotic dynamics and its resonance spectrum is well described by a GOE 



261 ] . i.e. using Eqs. 161 and [T01 with m = 1. We here consider the A = 0.125,0.150 billiards 



that exhibit mixed regular-chaotic dynamics, which is predominantly regular. We have 
performed a least-square analysis of the NNSD and LNV for these billiards using Eqs. [6] 
and [10J respectively, taking m as a real parameter. The best fit values are m = 3.21,2.62 
for the billiards with A = 0.125, 0.150, respectively. Figure 1 shows the result of comparison 
of the experimental NNSD for these billiards with Eq. El while the result for the LNV is 
given in FIG 2. We note,however, the interpretation of m as the number of spectra that 
are being superimposed suggests that it should be an integer. For this reason, we show in 
the two figures the results of calculation with 3 and 4 sequences for the A = 0.125 billiard 
and 2 and 3 sequences for the A = 0.150 one. The figures show that the agreement with the 
fractional value of m is not so much better than with the integer values of the parameter. 
In both cases, the parameter m can be taken equal to 3 for both the NNSD and LNV in 
spite of the fact that the NNSD is close to a Poisson distribution while the LNV shows a 
large amount of spectral rigidity. This unusual situation is in favor of the validity of the 
present model. To demonstrate this we show an analysis of the NNSD in FIG 3 and of the 
LNV in FIG 4 using the Berry- Robnik model 13j. The best-fit value of the parameter q 
that measures the fractional volume of the regular part of the phase space is found for the 
NNSD to be 0.585. This quite different from the value 0.156 that fits the LNV. We note 
that the agreement between the prediction of the Berry- Robnik model and the experimental 
LNV is worse than the agreement with our model. Concerning the NNSD, both models fail 
to describe the depletion in the number of events in the first bin. There is 100 spacings in 
this bin, so that the statistical error is 10 %. Thus the depletion is statistically significant. 
The disagreement reflects the partial neglect of level repulsion in both model where the 
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superimposed sequences are considered as independent. 

The expression for NNSD of a spectrum composed of independent sequences with non- 
equal fractional densities /j is more complicated. It has been shown in H] that the 

m 

NNSD in this case essentially depends on a single parameter, namely (/) = /?, which is 

i=i 

the mean fractional level density for the superimposed sequences; the statistical weight of 
each sequence is given again by its fractional density. For a superposition of equal sequences, 
fi = (/) = l/ m - Therefore, Eq. [6] can approximately be used to describe a superposition of 
independent but not equal sequences by considering m as a parameter, not necessarily taking 
integer values. The non-integer parameter m will play the role of an effective number of the 
constituting sequences m = 1/ (/) . One can adopt this interpretation of the parameter m 
if one sees that the fit in FIGs 1 and 2 are deteriorated by taking m = 3 instead of 3.21 or 
2.62, 

To summarize, we consider a model for systems with regular-chaotic dynamics in which 
the energy spectrum is represented by an independent sequences of levels, each one modeled 
by a Gaussian random ensemble. By varying the effective number of sequences, the model 
interpolates between the Poissonian spectrum for the regular system where the spectrum 
consists of infinite number of sequences and that of a chaotic system whose spacing distribu- 
tion is approximated by the Wigner surmise. We show that the model successfully describe 
both the NNSD and LNV for a nearly integrable Limagon billiard with the same value of 
the model parameter. 
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Figure Caption 

FIG 1. Comparison between the NNSD for two nearly integrable Limagon billiards with 
the same distribution for m independent GOE sequences. 

FIG 2. Comparison between the LNV for two nearly integrable Limagon billiard with 
the same distribution for m independent GOE sequences. 

FIG 3. Comparison between the NNSD for a nearly integrable Limagon billiard (A = 
0.125) with the same distribution calculated using the Berry-Robnik semiclassical method. 

FIG 4. Comparison between the LNV for a nearly integrable Limagon billiard (A = 0.125) 
with the same distribution calculated using the Berry-Robnik semiclassical method. 



10 




r 




r 



